Universality of efficiency at maximum power 
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We investigate the efficiency of power generation by thermo-chemical engines. For strong coupling 
between the particle and heat flows and in the presence of a left-right symmetry in the system, we 
demonstrate that the efficiency at maximum power displays universality up to quadratic order in 
the deviation from equilibrium. A maser model is presented to illustrate our argument. 
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The concept of Carnot efficiency is a cornerstone of 
thermodynamics. It states that the efficiency of a cyclic 
("Carnot") thermal engine that transforms an amount 
Q r of energy extracted from a heat reservoir at tem- 
perature T r into an amount of work W is at most 
77 = W/Q r < rj c = 1 — Ti/T r , where Tj is the temperature 
of a second, colder reservoir. The theoretical implications 
of this result are momentous, as they lie at the basis of the 
introduction by Clausius of the entropy as a state func- 
tion. The practical implications are more limited, since 
the upper limit r\ c ("Carnot efficiency") is only reached 
for engines that operate reversibly. As a result, when the 
efficiency is maximal, the output power is zero. By op- 
timizing the Carnot cycle with respect to power rather 
than efficiency, Curzon and Ahlborn found that the cor- 
responding efficiency is given by t\ca = 1 — y/Ti/T r [lj. 
They obtained this result for a specific model, using 
in addition the so-called endo-reversible approximation 
(i.e., neglecting the dissipation in the auxiliary, work pro- 
ducing entity). Subsequently, the validity of this result 
as an upper bound, as well as its universal character, 
were the subject of a longstanding debate. In the regime 
of linear response, more precisely to linear order in r) c , 
it was proven that the efficiency at maximum power is 
indeed limited by the Curzon-Ahlborn efficiency, which 
in this regime is exactly half of the Carnot efficiency, 
V < "Oca — Vc/2 + 0(r%) 0- The upper limit is reached 
for a specific class of models, namely, those for which the 
heat flux is strongly coupled to the work-generating flux. 
Interestingly, such strong coupling is also a prerequisite 
for open systems to achieve Carnot efficiency 0, |j| . In 
the nonlinear regime, no general result is known. Efficien- 
cies at maximum power, not only below but also above 
Curzon-Ahlborn efficiency, have been reported [!, H, 0, H] • 
However, it was also found, again in several strong cou- 
pling models 0, B that the efficiency at maximum 
power agrees with tjca up to quadratic order in r\ c , i.e., 
t] = t] c /2+t]^/8 + 0(7]^), again raising the question of uni- 
versality at least to this order. In this letter we prove that 
the coefficient 1/8 is indeed universal for strong coupling 



models that possess a left-right symmetry. Such a univer- 
sality is remarkable in view of the fact that most explicit 
macroscopic relationships, for example the symmetry of 
Onsagcr coefficients, arc limited to the regime of linear 
response. The interest in strong coupling is further moti- 
vated by the observation that it can naturally be achieved 
in nano-devices 2, II], [12] • To complement our theoret- 
ical discussion, we also present a detailed study of a ther- 
mal nano-machine based on the operation of a maser [l3| . 
It can be solved analytically and illustrates all the above 
mentioned features. Depending on the value of the Ein- 
stein coefficients, the efficiency of the maser at maximum 
power may be above or below Curzon-Ahlborn. However, 
when the Einstein coefficients are equal, the predicted 
universality is observed, with the universal value 1/2 for 
the linear coefficient, and the quadratic coefficient equal 
to 1/8. 

In view of the interest of our analysis for small scale 
systems, and in order to establish the connection with the 
subsequent discussion of the maser model, we derive the 
main results on the basis of a stochastic thermodynamic 
analysis as formulated for a master equation description 
of a driven open system [II], EH • As we will show in pass- 
ing, this formalism is fully consistent with macroscopic 
thermodynamics. 

The system under consideration is characterized by a 
set of states i of energy q and number of particles Ni . It 
exchanges particles and energy with two reservoirs fi = 
I, r, with inverse temperatures /3 M and chemical potentials 
fifj,, respectively. The probability of finding the system 
in state i at time t is denoted by pi (t) . The state of the 
system evolves in time according to a stochastic process 
which is described by the master equation 



(i) 



As a result of conservation of total probability, the 
stochastic rate matrix obeys the usual condition 
J2i Wij = 0. Wij is the probability per unit time to make 
a transition to state i from state j. We assume that these 
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transition rates are expressed as sums of independent 
contributions from the two reservoirs, Wij = 5Z„ 
{y = I, r). To reproduce the correct properties at equi- 
librium, it follows that each of the separate rate matrices 
satisfies detailed balance with respect to the grand 
canonical distribution at the prevailing temperature and 
chemical potential, 
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The average energy and matter currents entering the sys- 
tem from the reservoir v are given by 



l [ E\t)=Y,W% ) p j {t){e i -e j ), 



(3) 



(4) 



The rate of change of the total energy of the system and 
the (chemical) work per unit time on the system read 



£{t) = Y J Mt)z l = Y, I E ) ( t )' 
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The corresponding total average heat flow follows from 
energy conservation, Q(t) = Y, v Q^Kt) = £(*) - M*). 
In particular, the (average) heat flow from the reservoir 
v into the system is given by 



Q<r)(t) =1% ) 
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(6) 



The entropy of the system is taken to be the usual sys- 
tem entropy S(t) = — ^2 m p m (t)^Pm{t) (Boltzmann's 
constant ks — 1)- Using the master equation ([T]), one 
easily verifies that the rate of change of this entropy can 
be written in the form of a balance equation, namely, 
S(t) = Si(t) + S e (t). Here, 5»(t) is the non-negative total 
entropy production for the physical processes represented 
by the master equation, 



Si(t) 
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> 0. 



(7) 



Using Eq. ©, one verifies that the entropy flow into the 
system is given by the familiar thermodynamic expres- 
sion in terms of the heat fluxes, S e (t) = J2 V Q {u) {t)/T v . 

We focus on the case of a nonequilibrium steady state. 
From S = we have that Si — —S e . Also, current conser- 
vation at steady-state implies that ^2 V 1^' = 1 



M 



0. As a result, the entropy production can now be written 
in the traditional bilinear force-flux form 11611 



Si 



TmI^ > 0, 



(8) 



with the standard expressions for the thermodynamic 
forces, 
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We are interested in the operation of the device as a 
heat engine that carries particles uphill in the chemical 
potential, driven by a heat current from the hot to the 
cold reservoir. With no loss of generality, we henceforth 
assume that T r > T; and \x r < /i;. As mentioned before, 
we focus on the power generated by the device, which, 
due to current conservation in the steady-state, reads: 



v = -W = -(nr -m)f 



(r) 
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(10) 



The resulting efficiency of producing chemical work from 
the heat pumped out of the hot reservoir r is [see Eq. ([BJ] 
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The above formalism can be further simplified for the 
case of strong coupling between the energy and matter 
flux, defined as 



T 



el 
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(12) 



This condition implies that energy is exclusively trans- 
ported by particles of a given energy e. Such a selec- 
tion is quite natural in quantum nano-devices such as 
the maser 2 (see also below), and in thermo-electrical 
nano-devices [11|, but it can also occur in the classical 
context, for example for Kramers' escape, where the par- 
ticles that cross the barrier have precisely the minimum 
energy needed to do so [§]. Using Eq. (TTJ]), one can now 
rewrite the entropy production ([8]) in the simple form 
Si = Tl. Here 1 is the current introduced in Eq. (jX2j) 
and the associated thermodynamic force T can be ex- 
pressed in terms of dimensionless scaled energies xi and 



T = xi 
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XI 
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(13) 



Hence the two flows and forces collapse into a single flux 
1 and a single corresponding thermodynamic force T ', 
respectively. Note that equilibrium, that is, zero entropy 
production, is reached for T = 0. This does not require 
that the forces Te and Tm be zero separately. In fact in 
the vicinity of T = the device can operate at Carnot 
efficiency, see for example 0, li, 11 1. Using Eq. (Tl2|l . the 
power PH| becomes 

V = -W= {T r x r ~ Tixi)l = -T r {T - rj c xi)l, (14) 

and the thermodynamic efficiency (jTTJ) reads 

W - /^ , \ x l ,,,, 

T)=— =l-(l-7?c) — . (15) 
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The properties of the system are contained in the de- 
pendence of the flux X on the variables xi and x r , 
T = T(x r ,Xi). Its explicit expression is obtained by 
inserting the steady state solution of the master equa- 
tion |T]) into the expression ^ for the mass flux. 

To identify the regime of maximum power, we proceed 
in two steps. The extremum of power with respect to J- 
is determined by the condition 



T. 



I + - x lVc )drl = Q. 



(16) 



Since we are interested in the behavior around equilib- 
rium, including the first nonlinear correction term to the 
linear regime, we expand the current T{x ri x{) — T(xi — 
xi) to quadratic order in T, 1 = LT + MT 2 + 0(F 3 ). 
The Onsager coefficients are given by L — —2[(xi,xi) 
and M = I'^xi, xi)/2. The primes denote the num- 
ber of derivatives and the sub-indices indicate whether 
these derivatives are taken with respect to the first or sec- 
ond variable. Furthermore, since T has to become zero 
when r\ c goes to zero, we can write, again to quadratic 
order, that T = b\r\ c + cirj 2 + 0(rfc). Insertion in the 
extremum condition (|16j) allows us to identify the coef- 
ficients b\ — xi/2 and c\ — Mxf/(8L). The resulting 
expression for the efficiency (1151) reads 



V 



V. 



1 Mxi 



2 v 4 8L 
Next, we maximize power with respect to xi, 



dxi 



-r) c l + [T - xi7] c )d xi l = 0. 



(17) 



(18) 



It suffices to find the result to lowest order in rj c . Inserting 
the expansions 1 = LT+O^ 2 ) and T = r] c xi/2 + 0(rj 2 ), 
one finds that xi — —2L/d Xl L. Combined with Eq. (TTT)) . 
we finally arrive at the following result for the efficiency 
at maximum power, valid up to quadratic order in rj c : 



M 



L' 4 



(19) 



We conclude that, while we recover the universal value of 
the coefficient 1/2 in the linear term, the coefficient of the 
quadratic term is in general model dependent. However, 
as we now proceed to show, the appearance of the coeffi- 
cient 1/8 in several previously studied models 0,11,0] de- 
rives from the fact that these models possess a left-right 
symmetry. More precisely, such a symmetry implies that 
the switching of temperatures /3„ and chemical potentials 
[i v leads to an inversion of the flux, 



T(x r ,Xl) = —J(xi,X r ). 



(20) 



By deriving both sides with respect to x r and xi and then 
setting x r — xi, one finds that 2" 2 {xi : xi) = 0. Together 
with d Xl L = —T'l-yixi^xi) — T'^ixi, xi), we conclude that 



FIG. 1: Illustration of the maser model 

the condition 2M = —d Xl L is verified and universality 
of the coefficient 1/8 is established under the symmetry 
specified in Eq. PO]) . 

To illustrate these findings, we turn to the analysis 
of the maser model introduced in 13j, see Fig. [TJ The 
system possesses three energy levels e.;, i = 1,2,3. It ex- 
changes photons with three equilibrium black bodies R, 
L and S (temperatures T r , TJ and T s ) with correspond- 
ing specific frequencies hv r = 63—61, hv\ = £3 — £2 and 
hv s = 62 — ei- The reservoirs R, L and S control the tran- 
sitions 1 — 3, 2 — 3 and 1 — 2, respectively. The stochastic 
dynamics of these transitions are described by the master 
equation ([1]) , with rates corresponding to the processes of 
absorption, spontaneous emission, and stimulated emis- 
sion of the photons. They are given by — T r n(x r ) 
(absorption of a photon from R) and = T r [l+n(x r )] 
(spontaneous and stimulated emission of a photon into 
R), and identical expressions for transitions 2 — 3 and 
1 — 2, with the indices r replaced by I and s respec- 
tively. Here, we introduced the Bose-Einstein distribu- 
tion n(x) = [cxp(x) — 1] , with the scaled energies 
x r = hv r /T r , xi = hvi/Ti and x s = hv s /T s , and the 
reduced Einstein coefficients T v . 

To transform the system into a thermal engine, we 
consider the high temperature limit T s — » 00 (x s — > 0). 
The reservoir S effectively becomes a repository of work, 
since heat stored in a reservoir at infinite temperature 
can be recuperated at 100% (the corresponding Carnot 
efficiency being equal to 1). We next note the cyclic na- 
ture of the transitions: Starting from state 1, in order 
to deposit the amount — w — hv s as work into the S- 
reservoir (transition 2 — * 1), the system first needs to 
absorb a photon q( r > = hv r from the hot reservoir (tran- 



sition 1 — > 3) and next deposit —q 



hvi into the 



cold reservoir (transition 3 — > 2). This (with the reverse 
process) is the only available cycle. The corresponding 
efficiency of the cycle reads: 



V 



-w 



= !-(!-%) — ■ (21) 

e 3 - ei x r 



At steady state, the system will, on average, run through 
X such cycles per unit time, with corresponding heat flows 
Qy 1 — q r J and Q} 1 ' = q{I. The power of the device is 
given by 



-wT - 



)(r) 



Q 



3(0 



(T r a 



TixAl. 



(22) 
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FIG. 2: (Color online) (A) Efficiency at maximum power 
compared with Carnot efficiency (straight dashed line) and 
Curzon-Ahlborn efficiency (dotted line). (B) Scaled energies 
Xi and x r and maximum power V (note that xi and x r only 
depend on the ratio I^/IV but V not). 



We thus recover the previously derived results of the 
strong coupling regime, cf. Eqs. (fT4"|) and (fTS)) . To com- 
plete the analysis, we need to evaluate the steady-state 
current T. This is a matter of algebra, involving the 
steady state solution of the master equation using the 
transition rates given above. One finds: 



1 = 



(e Xi 



(1 + 2e x ')(e x - - l)r z + (1 + 2e x -)(e x ' - l)T r 



•(23) 



Concerning efficiency of the device at maximum power, 
we can now invoke the general conclusions mentioned 
earlier. The symmetry criterion (|20[) for the current is 
only satisfied when T; = T r . Under this condition, the 
efficiency at maximum power displays the universal co- 
efficient 1/8 for the quadratic term, in addition to the 
universal linear coefficient 1/2. This observation is con- 
firmed by an explicit calculation for the model under con- 
sideration. We find: 



v = ■£ + 1 



3(r ; - r r 



(F ; +r r )(3cosha-|-sinhcv)y' 8 



where a = 1.77676, solution of a transcendental equation 
2 + a + 2e a + 2(a — 2)e 2a — 0, is also the asymptotic value 
of xi and x r when r\ c — > 0. To complete the picture, we 



have reproduced, in Fig.[2jA), the efficiency at maximum 
power as a function of rj, with r\ £ [0, 1], for the cases 
r;/r r = 0,1 and oo [l3]- All three curves are remark- 
ably close to the Curzon-Ahlborn efficiency, even though 
the efficiency is slightly larger in the first two cases and 
slightly less in the last case. In view of their technological 
interest, we also include in Fig. [2fB) the corresponding 
maximum power and the operational conditions of the 
scaled energies xi and x r . 
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